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We  report  the  observation  of  high  resolution  (0.008  cm'*)  infrared  absorption  spectra  of 
CO2  molecules  isolated  in  solid  parahydrogen  (PH2)  matrices  at  T  =  2.4  to  4.8  K.  Several 
extremely  sharp  (0.01  to  0.04  cm'*  fiill-width-at-half-maximum)  absorption  features  appear 
in  the  2343.5  to  2345  cm'*  region.  We  assign  the  three  strongest  peaks  to  the  V3  mode  of 
isolated  CO2  molecules.  The  spectra  are  consistent  with  trapping  of  the  CO2  molecules  in 
three  distinct  double-substitutional  sites  in  hep  and  fee  regions  of  the  pH2  solid.  We  offer 
several  hypotheses  as  to  the  origins  of  the  numerous  weaker  absorption  features. 


I.  Introduction 

This  manuscript  is  a  status  report  on  our 
continuing  effort  to  observe  and  analyze  vibrational 
spectra  of  molecular  dopants  isolated  in  cryogenic  solid 
parahydrogen  (PH2)  hosts.  As  detailed  below,  our 
investigation  of  the  microscopic  structures  and 
dynamics  underlying  many  of  these  spectra  is  at  a 
preliminary  stage.  Despite  this,  we  believe  the 
experimental  observations  themselves  merit  being 
reported  at  this  time.  We  hope  that  these  new  high- 
resolution  data  will  stimulate  interest  in  the  rigorous 
testing  of  theoretical  models  of  vibrational 
spectroscopy  of  impurities  in  condensed  phases. 

The  suitability  of  solid  pH2  as  a  host  for  high- 
resolution  (Av/v  ~  10'^)  vibrational  spectroscopy  was 
discovered  by  Oka  and  co-workers  [1-5],  and  further 
investigated  by  Shida,  Momose  and  co-workers  [6-10], 
and  by  Winnewisser  and  co-workers  [1 1-14].  The 
favorable  properties  of  solid  pH2  contributing  to  this 
phenomenon  have  been  discussed  in  detail  [2,3,9]. 

Here  we  note  simply  that  solid  pH2  is  very  "forgiving" 
of  highly  non-equilibrium  sample  preparation 
techniques  such  as  direct  gas-to-solid  condensation. 
Structural  defects  in  such  samples  may  be  less 
numerous,  and/or  have  a  smaller  influence  on  dopant 
vibrational  spectra,  than  in  other  vapor  deposited 
cryogenic  van  der  Waals  solids  traditionally  employed 
as  hosts  for  matrix  isolation  spectroscopy  (^S). 

Our  work  on  solid  pH2  is  supportedoy  the  U.S. 
Air  Force's  High  Energy  Density  Matter  (HEDM) 
program  [15];  our  project's  ultimate  objective  is  to 
demonstrate  practical  energy  storage  in  cryogenic 


solids  for  use  as  advanced  chemical  rocket  propellants. 
Thus,  our  attention  is  focused  on  the  problems  of  (a) 
trapping  large  concentrations  of  isolated  energetic 
dopants,  (b)  scaling-up  sample  quantities  beyond  the 
thin  films  typically  encountered  in  MIS  studies,  (c) 
rigorously  characterizing  the  chemical  identities, 
concentrations,  and  trapping  site  structures  of  the 
incorporated  energetic  species,  and  (d)  evaluating  the 
thermal  and  chemical  stabilities  of  these  prototypical 
"cryosolid"  propellants.  These  last  two  tasks  motivate 
our  interest  in  the  spectroscopy  of  dopants  in  solid  pH2. 

Our  new  "rapid  vapor  deposition"  technique 
produces  easily  doped,  millimeters-thick  pH2  solids  of 
remarkable  optical  clarity  [16,17].  The  compatibility  of 
this  method  with  most  MIS  dopant  trapping  schemes 
enables  the  isolation  in  solid  pH2  of  a  wide  variety  of 
previously  inaccessible  chemical  species.  The  excellent 
transparency  of  rapid  vapor  deposited  samples  is 
unusual  for  pH2  solids  produced  by  vapor  deposition 
onto  a  substrate-in- vacuum  [18],  and  enables  their 
characterization  by  optical  methods. 

Our  early  work  on  doped  cryogenic  hydrogen 
solids  employed  medium-resolution  ultraviolet/visible 
and  infrared  (IR)  diagnostics  [19-21].  However,  we 
recently  added  a  high-resolution  IR  capability  and  were 
repaid  by  the  observation  of  very  sharp  (~  0.01  cm  * 
full-width-at-half-maximum,  FWHM)  vibrational 
absorption  features  for  dopants  in  rapid  vapor  deposited 
samples.  Moreover,  high-resolution  IR  spectra  of 
several  (presumably)  simple  molecular  dopants  in  solid 
pH2  exhibit  amazingly  rich  structure.  Not  only  are  the 
vibrational  transition  energies  often  more  precisely 
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defined  than  for  the  same  dopants  in  rare  gas  solid 
(RGS)  hosts,  but  such  spectra  also  typically  include  a 
larger  number  of  distinct,  well  resolved  features.  The 
increased  precision  and  complexity  of  these  spectral 
measurements  may  require  extraordinary  efforts  for 
their  complete  assignment,  but  the  reward  should  be  a 
deeper  and  more  quantitative  understanding  of  dopant 
trapping  site  structures  and  vibrational  dynamics  in 
solid  pH2  than  has  been  achieved  to  date  by  MIS  in 
RGS  hosts. 

We  are  in  the  process  of  surveying  high- 
resolution  IR  absorption  spectra  for  a  variety  of  dopants 
in  solid  pHi.  A  substantial  database  of  such  spectra 
should  help  to  determine  which  physical  effects 
dominate  the  vibrational  spectroscopy  of  different 
classes  of  dopants  in  solid  pH2.  For  example,  spectra  of 
CH4/PH2  solids  [6,7,9,10,22],  a  case  in  which  the 
spherical-top  dopant  fits  easily  into  a  single- 
substitutional  trapping  site,  accordingly  show  line 
spacings  and  intensity  progressions  commensurate  with 
slightly  hindered  rotors  [6,22,23-27].  Structurally,  the 
spectra  demonstrate  that  in  as-prepared  rapid  vapor 
deposited  pH2  solids  some  of  the  CH,  molecules  are 
initially  trapped  in  face-centered-cubic  (fee)  regions, 
and  that  these  regions  convert  irreversibly  to 
hexagonal-close-packed  (hep)  upon  annealing  [22].  The 
polarization  dependences  of  the  CH4/PH2  transitions 
further  indicate  that  annealing  also  results  in  improved 
alignment  of  the  hep  crystallites'  c-axis  orientations 
with  the  deposition  substrate  surface  normal  [22]. 
Dynamically,  all  the  CH4/PH2  spectra  can  be  rigorously 
assigned  by  considering  the  rotational-vibrational 
dynamics  of  CH4  at  the  center  of  static  external  fields 
having  Dah  [6]  and  Oh  [22]  symmetries.  Thus,  a 
satisfactory  explanation  of  these  transition  energies 
does  not  require  consideration  of  constrained  rotational- 
translational-coupling  (RTC)  in  a  rigid  trapping  site 
[28-31],  or  of  more  general  RTC  interactions  within  a 
dynamic  trapping  site  structure  [32-37].  These  lessons 
learned  from  the  CH4/PH2  studies  provide  an  excellent 
starting  point  for  the  assignment  of  our  recently 
observed  highly  structured  CO/pH2  spectra  [38], 
another  case  in  which  the  dopant  molecule  should  fit 
readily  into  a  single  substitutional  vacancy. 

In  stark  contrast,  we  find  that  spectra  of  many 
larger  species,  such  as  lower  symmetry  polyatomic 
molecules,  and  hydrogen-  or  van  der  Waals-bonded 
clusters  of  smaller  dopant  molecules,  show  no  evidence 
of  overall  rotation.  We  believe  end-over-end  rotation  is 
inhibited  by  strong  anisotropic  interactions  within  the 
lower  symmetry  multi-substitutional  trapping  sites 
required  to  accommodate  such  large  species.  Thus, 
these  larger  species  likely  exist  in  the  "librational" 
limit,  oscillating  around  certain  equilibrium  orientations 


instead  of  rotating  as  a  unit  [24].  Yet,  despite  these 
strong  dopant-host  interactions,  many  of  these  larger 
dopant  systems  exhibit  even  sharper  vibrational 
absorptions  than  those  observed  for  smaller  dopants 
trapped  in  highly  symmetrical  single  substitutional 
sites!  A  satisfactory  analysis  of  the  spectra  of  these 
larger  librating  systems  must  deal  with  the  additional 
complication  of  simultaneously  considering  both 
inherent  and  trapping-site-induced  spectral  features. 

Before  attempting  such  analyses,  we  turn  first 
to  a  simpler  model  system  of  a  linear  triatomic 
molecule  such  as  CO2  in  solid  pH2.  Our  IR  spectra  of 
CO2/PH2  samples  show  complicated  patterns  of  sharp 
lines  which  do  not  appear  to  match  a  rotational 
dynamics  {vide  infra),  and  so  may  be  comparable  in 
this  respect  to  the  larger  dopants.  However,  the  gas- 
phase  vibrational  transitions  of  CO2  are  very  well 
known;  thus,  any  peculiar  features  in  spectra  of 
CO2/PH2  are  immediately  attributable  to  either  dopant- 
dopant  interactions  or  to  interactions  with  the  pH2  host. 
We  think  the  CO2/PH2  system  will  prove  to  be  an 
excellent  arena  for  evaluating  the  various  elements  of 
the  RTC  models  in  the  limit  of  strongly  hindered 
rotation. 

In  what  follows  we  describe  briefly  our 
experimental  methods  and  present  our  preliminary 
high-resolution  IR  absorption  data  for  the  CO2/PH2 
system.  We  will  show  that  the  strongest  absorption 
features  are  due  to  isolated  CO2  molecules,  and  not  to 
(C02)n  clusters.  We  discuss  these  monomer  spectra  in 
terms  of  candidate  trapping  site  structures,  and  the 
corresponding  rotational  vs.  librational  dynamics  of  the 
trapped  molecules.  We  will  give  a  list  of  hypotheses 
about  the  possible  origins  of  the  numerous  weaker 
absorption  features,  but  attempt  no  quantitative 
application  of  crystal  field  or  RTC  models  at  this  time. 
We  end  this  manuscript  with  a  summary,  and  list  some 
of  our  immediate  plans  for  future  experimental  and 
theoreticarwork  to  test  the  proposed  hypotheses. 

n.  Experimental 

Our  experimental  apparatus  and  sample 
preparation  tectoiques  have  been  described  in  detail 
before  [16,17,21].  Doped  pH2  solids  are  prepared  by 
rapid  vapor  deposition  of  precooled  pH2  gas  and  room 
temperature  CO2  gas  onto  a  BaF2  substrate  cooled  to  T 
»  2  K  in  a  liquid  helium  (IHe)  bath  cryostat.  We 
operate  the  ortho/para  hydrogen  converter  at  15  K, 
yielding  a  flow  of  pre-cooled  pH2  containing  «  0.01  % 
residual  orthohydrogen  (0H2)  impurities  [17].  The  pH2 
flow  impinges  upon  the  substrate  at  a  45°  angle;  during 
a  deposition  the  pressure  of  uncondensed  pH2  gas 
remains  below  ~  1 0''*  Torr.  The  CO2  dopant  is  metered 
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Figure  1.  IR  absorption  spectra  of  three  as-deposited  pH2 
solids  at  T  =  2.4  K  containing  different  concentrations  of 
natural  isotopic  abundance  CO2.  Trace  (la)  is  for  a  13  PPM 
CO/pH2  sample  containing  ~  0,01  PPM  CO2  as  an 
unintentional  impurity;  sample  thickness  is  2.7  mm.  Trace 
(lb)  is  for  an  8  PPM  HCl/pH2  sample  containing  «  0.04  PPM 
CO2  as  an  unintentional  impurity;  sample  thickness  is  3.0 
mm.  Trace  (Ic)  is  for  a  1.2  PPM  CO2/PH2  sample  that  is  1.4 
mm  thick.  Hie  inverted  triangles  at  the  top  of  die  figure 
indicate  the  positions  of  several  gas  phase  ro-vibrational  lines 
for  the  V3  mode  of  *^C^^02  [40].  Trace  (la)  has  been  rescaled 
by  a  multiplicative  factor  of  4;  all  the  traces  have  been 
displaced  vertically  for  ease  of  presentation. 

into  the  deposition  cryostat  from  a  separate  gas 
handling  manifold;  the  dopant  inlet  is  situated  at  an 
angle  of  45^  from  the  substrate  surface  normal,  at  90° 
to  the  pH2  source.  Individual  sample  preparation 
details  are  given  belo>v  in  the  figure  captions. 

We  record  IR  absorption  spectra  of  our 
CO2/PH2  samples  across  the  800  to  7800  cm’’  range  at  a 
resolution  of  0.008  cm'*;  the  main  optical  axis  is 
parallel  to  the  substrate  normal.  The  Fourier  Transform 
IR  spectrometer  (Bruker  IFS120HR)  is  equipped  with  a 
glowbar  source,  a  KBr  beamsplitter,  and  a  liquid 
nitrogen  cooled  HgCdTe  detector.  To  accommodate 
the  ER  diagnostic,  the  entire  optical  path  is  enclosed 
within  a  0.5  m^  polycarbonate  box  purged  with  a 
constant  flow  of  dry  N2  gas. 


I  ■  ■  ■  ■  I 

P 


2343.0  2343.5  2344.0  2344.5  2345.0  2345.5  2346.0 

wavenumber  (cm'h 

Figure  2.  Annealing  behavior  of  the  1,2  PPM  CO2/PH2 
sample  described  in  figure  1.  Trace  (2a)  is  an  expanded  view 
of  the  as-deposited  spectrum  presented  above  in  trace  (Ic). 
The  temperature  sequence  is:  (2a)  T  =  2.4  K,  (2b)  T  =  4.8  K, 
(2c)  T  =  2.4  K,  (2d)  T  =  4.8  K,  and  (2e)  T  -  2.4  K.  The 
labels  a'  and  a”  denote  two  new  features  that  grow  in  during 
annealing. 

m.  Results 

Figure  1  shows  the  V3  region  of 

absorption  spectra  from  three  different  as-deposited  pH2 
solids  containing  CO2  concentrations  ranging  from 
0,01  PPM  up  to  1.2  PPM.  The  concentrations  are 
estimated  as  described  previously  [21]  using  a  value  of 
550  km/mol  for  the  V3  CO2  integrated  absorption 
coefficient  [39].  The  peaks  labeled  a,  p,  and  y 
dominate  all  three  spectra,  but  show  minor  sample-to- 
sample  variations  in  relative  intensities.  Not  shown  are 
the  V3  ^^C^^02  and  V3  regions  of  the  spectrum 

depicted  in  trace  (Ic),  in  which  similar  features,  also 
matrix- shifted  by  «  -5  cm'^  from  the  gas-phase 
vibrational  band  origins  [40,41],  appear  for  each  natural 
abundance  isotopomer.  The  observed  peak  positions 
and  linewidths  are  summarized  in  Table  I. 
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Table  1.  Peak  positions  (cm  *)  and  widths  (FWHM,  rounded  to  nearest  0.005  cm  ),  and  assignments  for  the  main  IR  ^ 
absorptions  in  CO2/PH2.  The  labels  a,  p,  and  y  refer  to  the  three  largest  peaks  observed  in  as-deposited  samples;  the  labels  a' 
and  a"  refer  to  two  additional  strong  features  that  appear  between  a  and  p  upon  temperature  cycling.  The  "as-deposited" 
spectra  are  taken  at  T  =  2.4  K,  the  "annealing"  spectra  at  T  =  4.8  K,  and  the  "annealed"  spectra  upon  re-coolmg  to  T  -  2.4  K. 
The  gas-phase  vibrational  band  origins,  vo,  are  estimated  from  data  in  reference  40. 


- peak  positions  (widths) - 

-  - 

label 

as-deoosited 

annealine 

annealed 

assignment 

a 

2278.296  (0.040) 

2278.300(0.015) 

2278.299  (0.010) 

V3  ‘**€**02 

a' 

n.o. 

Vo  =  2283 .48 

a” 

2278.454W  (0.025) 

p 

2278.861  (0.030) 

2278.872  (0.030) 

2278.867  (0.025) 

Y 

2279.099  (0.015) 

a 

2326.835  (0.035) 

2326.842(0.015) 

2326.842(0.015) 

V3  •«o*'c*’’o 

a' 

n.o. 

Vo  =  2332.11 

a" 

n.o. 

6 

2327.4  Iw 

2327.42W 

2327.399  (0.010) 

r 

2327.422(0.010) 

Y 

2327.656  (0.020) 

a 

2343.813  (0.040) 

2343.818(0.020) 

2343.814(0.015) 

V3  ‘^C*®02 

a’ 

2343.912(0.015) 

Vo  =  2349.14 

a" 

2343.979(0.020) 

P 

2344.395  (0.025) 

2344.404  (0.030) 

2344.397  (0.040)* 

Y 

2344.639  (0.020) 

w  -  weak,  signal/noise  ~  1 

n.o. not  observed,  signal/noise  <  1 

*  -  shoulder  at  2344.372  cm'* 


Figure  2  shows  the  effects  of  repeated 
temperature  cycling  on  the  V3  *^C**02  band  of  the  1.2 
PPM  CO2/PH2  sample  depicted  in  trace  (Ic).  The  most 
pronounced  change  observed  during  the  initial 
wanning  from  T  =  2.4  to  4.8  K  is  the  strong 
irreversible  decrease  in  the  intensity  of  the  y  peak  at 
2344.64  cm'*;  the  same  behavior  is  observed  for  the 
other  CO2  isotopomers,  as  well.  We  also  note  the 
simultaneous  growth  of  the  a'  peak  at  2343.91  cm' , 
the  weakening  and  sharpening  of  the  a  line,  the 
strengthening  and  broadening  of  the  P  peak,  and  the 
appearance  of  a  »  0.5  cm  *  broad  absorption  lump  near 
2344.6  cm'*.  Trace  (2c)  shows  that  upon  re-cooling  to 
T  =  2.4  K  the  a  and  a'  peaks  weaken,  a  new  a"  line 
appears  at  2343.98  cm'*  together  with  other  weaker 
features  in  the  2343.6  to  2343.9  cm '  region,  the  P 
feature  strengthens  and  broadens  further,  and  the 
2344.6  cm'*  lump  broadens  and/or  diminisnes. 
Repeating  the  temperature  cycle,  traces  (2d)  and  (2e), 
shows  the  largely  reversible  nature  of  these  changes. 


Figure  3.  Continuation  of  the  spectra  presented  in  figure  2. 
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As  reported  in  Table  I,  the  broadening  of  the  V3 
*^C‘®02  P  feature  upon  annealing  is  largely  due  to  the 
appearance  of  a  shoulder  on  the  red  side  of  the  main 
peak;  the  V3  '*0‘^C'*0  P  feature  actually  splits 
completely  into  two  sharp,  well-resolved  peaks 
separated  by  «  0.02  cm'*. 

Figure  3  shows  the  2345  to  2350  cm'*  region 
of  the  same  spectra  as  depicted  in  figure  2.  Several 
sharp  new  peaks  appear  upon  temperature  cycling, 
most  notably  those  near  2345.25, 2345.37, 2346.67, 
and  2347.25  cm'*. 

IV.  Discussion 

A.  Isolation  of  CO2  Molecules 

We  assign  the  a,  p,  and  y  features  to  isolated 
CO2  molecules  as  opposed  to  (C02)n  clusters  due  to: 

(1)  the  extreme  dilutions  of  the  samples,  (2)  the  lack  of 
any  systematic  concentration  dependence,  and  (3)  the 
absence  (for  as-deposited  CO2/PH2  solids)  of  any  other 
absorption  features  outside  a  narrow  »  1  cm'*  window. 
Assessing  the  validity  of  this  last  point  requires  a  brief 
review  of  the  literature  on  (C02)n  clusters. 

The  gas-phase  vibrational  band  origin  for  the 
localized  asymmetric  stretch  mode  of  "slipped 
parallel"  CO2  dimers  is  observed  [42]  to  be  shifted  by 
+1.63  cm'*  relative  to  the  monomer  V3  *^C*®02  band 
origin  at  2349.14  cm'*;  theoretical  calculations  [43-45] 
predict  shifts  of  between  0  and  +2  cm'*,  depending  on 
the  dimer  geometry.  The  cyclic  Csh  symmetry  CO2 
trimer  band  origin  is  shifted  by  +2.58  cm  *  [46],  and 
the  nonrcyclic  C2  symmetry  CO2  trimer  shows  two 
bands  with  origins  shifted  by  -5.85  cm"  and 
+3.58  cm'*  [47],  all  shifts  again  referenced  to  the  gas- 
phase  monomer  V3  band  origin.  Larger  gas-phase 
(C02)n  clusters  [48]  and  ultra-fine  particles  [49]  show 
broad  absorptions  across  the  2340  to  2380  cm  *  range. 
The  absorption  maximum  in  pure,  crystalline,  natural 
isotopic  abundance  solid  CO2  occurs  at  2344.8  cm  , 
but  solid  samples  prepared  by  different  methods  can 
show  absorptions  anywhere  between  2330  and  2380 
cm'*  [50-53]. 

Matrix  isolation  studies  involving  general  X- 
CO2  complexes  are  far  too  numerous  to  list  here 
exhaustively  [54-56].  The  IR  absorption  spectrum  of  a 
1  PPM  CO2/N2  solid  carefully  prepared  from  the  melt 
shows  a  single  sharp  line  with  a  width  of »  0.007  cm 
at  T  =  1 1  K,  which  is  assigned  to  librating  CO2 
molecules  isolated  in  single-substitutional  sites  [57]. 
All  other  previously  reported  V3  region  spectra  of 
solely  CO2  doped  rare  gas  [58-62],  deutenum  [63,64], 
and  pH2  [65]  matrices  show  multiple  absorption 
features  spanning  a  5  to  10  cm  *  range.  These  features 
are  preparation,  concentration,  and  annealing 


dependent,  and  are  assigned  to  CO2  molecules  in 
multiple  trapping  sites,  and/or  to  aggregated  (C02)n 
species.  Most  directly  relevant  to  this  discussion  are 
the  spectra  of  (C02)n  clusters  in  solid  pH2  which  show 
a  complicated  pattern  of  peaks  throughout  the  2345  to 
2350  cm'*  region  [65]. 

Finally,  the  absence  of  additional  features  in 
traces  (la)  and  (lb),  for  which  the  CO2  is  actually  an 
unintentional  impurity  in  samples  containing  ~  10 
PPM  of  other  dopants,  further  supports  the  argument 
for  negligible  dopant  aggregation  under  the  present 
sample  preparation  conditions.  Not  shown  are  spectra 
from  other  experiments  at  higher  HCl  concentrations 
which  do  show  absorption  features  in  the  2346  to  2349 
cm'*  region  attributable  to  (HCl)ni(C02)n  complexes 
[66].  These  2  to  5  cm'*  blue  shifts  firom  the  matrix 
isolated  monomer  absorptions  are  in  line  with  the  +2 
cm'*  shift  reported  in  previous  HCl/C02/Ar  studies 
[67,68].  Thus,  clustering  during  deposition  to  form 
(C02)n  should  be  even  less  important  for  the  more 
dilute  1.2  PPM  CO2  sample  depicted  in  trace  (Ic). 

We  thus  conclude  with  confidence  that  the  a, 

P,  and  y  features  are  due  to  isolated  CO2  molecules. 

The  appearance  of  all  detectable  absorption  features  in 
as-deposited  samples  within  a  narrow  »  1  cm  *  region, 
especially  traces  (Ic)  and  (3  a),  further  suggests  that 
even  these  weaker  features  may  also  be  due  to  isolated 
CO2  molecules.  The  new  weak  features  which  appear 
only  upon  annealing,  especially  those  with  significant 
shifts  from  the  band  center,  may  be  due  to  (C02)n 
clusters.  We  will  discuss  these  possibilities  further  in 
the  following  sections. 

B.  CO2  Molecule  Trapping  Site(s) 

We  attribute  the  appearance  of  (at  least)  three 
separate  absorption  peaks  for  monomeric  CO2  to 
trapping  of  CO2  molecules  in  multiple  trapping  sites. 
Each  distinct  trapping  site  by  definition  corresponds  to 
a  different  time-averaged  structure  of  dopant  and  host 
molecules  and  so,  in  principle,  can  generate  a  distinct 
pattern  of  gas-to-matrix  spectral  shifts  for  the  dopant 
vibrational  transitions.  For  purposes  of  this 
discussion,  we  can  ascribe  the  formation  of  multiple 
trapping  sites  during  the  highly  non-equilibrium . 
sample  deposition  process  to  (1)  imperfections  in  the 
underlying  host  crystal  structure,  and/or  (2) 
mismatches  in  the  "sizes"  of  host  and  ^est  molecules. 

As  mentioned  in  the  Introduction,  the 
microscopic  structure  of  our  rapid  vapor  deposited  pH2 
solids  is  far  from  that  of  the  perfect,  single  hep  crystals 
that  can  be  produced  by  slowly  freezing  liquid  pH2 
[69].  The  hep  (...ABABAB...)  and  fee 
(...^CABCABC...)  structures  are  only  two  of  the 
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Figure  4.  Hypothetical  trapping  sites  for  CO2  molecules  in  solid  pH2  based  on  single-,  double-,  and  triple-substitutional 
vacancies  in  close-packed  planes  of  pH2  molecules. 


infinite  number  of  densest  close-packed  structures  that 
can  be  produced  by  stacking  of  close-packed  ("basal") 
planes.  Patterns  with  longer  repeat  units  (polytypism), 
non-periodic  structures  such  as  twin  and  stacking  faults, 
and  even  random  stacked  close-packed  structures  are 
also  possible  [70]. 

Fortunately,  the  combination  of  IR  and  Raman 
spectroscopies  of  pure  solid  pH2  reveals  the  symmetries 
of  the  sites  occupied  by  the  pH2  molecules  [71-74].  We 
have  not  detected  any  vacancy  defects  in  as-deposited 
undoped  pH2  solids,  and  have  further  shown  that  these 
samples  are  not  amorphous;  rather,  they  appear  to  be 
densest  close-packed  solids  made  up  of  separate  hep 
and  fee  regions  [16,22,74].  The  introduction  of  up  to  -s' 
1000  PPM  concentrations  of  dopants  like  CH4  and  CO, 
which  are  nearly  the  same  "size"  as  a  pH2  molecule  in 
the  van  der  Waals  sense,  does  not  significantly  affect 
the  pH2  structure.  Thus,  trapping  of  such  dopants  in 
single-substitutional  vacancies  results  in  the  formation 
of  two  distinct  classes  of  trapping  sites:  of  Dsh 
symmetry  in  hep  regions,  and  Oh  symmetry  in  fee 
regions.  The  experimental  CH4/PH2  and  CO/pH2 
spectra  are  all  consistent  with  this  relatively  simple 
picture  [22,38]. 

However,  if  our  spectroscopic  diagnostics  are 
only  sensitive  to  nearest  neighbor  PH2-PH2  interactions, 
then  we  would  necessarily  detect  only  two  pH2 
environments  for  any  arbitrary  stacking  pattern  of  basal 
planes,  i.e.\  sites  in  layers  (B)  embedded  in  regions  of 
local  hep  stacking  (XABAX,  “hep-like”),  or  local  fee 
stacking  (XABCX,  “fcc-like”).  Dopants  substituted 
into  hep-like  sites  would  experience  external  fields  of 
very  nearly  Dsh  symmetry,  those  in  fcc-like  sites  would 
exist  in  nearly  octahedral  environments.  Depending  on 
the  particulars  of  the  dopant-pH2  interactions,  the 


distinction  between,  for  example,  "hep-like"  and  "true 
hep"  may  be  undetectable  in  a  given  spectroscopic 
experiment. 

In  the  case  of  CO2  as  the  dopant,  the  situation  is 
further  complicated  by  the  larger  size  of  CO2  relative  to 
the  pH2  host.  Figure  4  shows  this  size  discrepancy  in  a 
semi-quantitative  manner,  with  the  CO2  molecule 
depicted  in  single-,  double-,  and  triple-substitutional 
vacancies  in  basal  planes  of  spherical  pH2  molecules. 
Overlap  between  molecular  outlines  indicates  repulsive 
interactions,  whereas  intervening  space  indicates 
attractive  interactions. 

The  pH2  molecules  are  drawn  as  circles  with 
diameters  representing  the  3.8  A  nearest-neighbor 
separation  in  solid  pH2  at  IHe  temperatures  [69].  We 
use  the  nearest  neighbor  separation,  instead  of  the 
minimum  separation  for  the  pH2-pH2  interaction 
potential,  to  include  the  effects  of  quantum  zero  point 
motion  on  the  structure  of  solid  pH2  which  would  be 
absent  in  a  simple  potential  energy  minimization.  We 
were  unable  to  find  an  angle-dependent  CO2-PH2 
interaction  potential  in  the  literature,  so  we  resorted  to  a 
recently  calculated  C02-He  potential  instead  [75]. 

Since  the  minima  of  the  (spherically  averaged)  He-He, 
H2-H2,  and  He-H2  potentials  occur  at  separations  of: 
2.94  A  [76],  3.4  A  [69],  and  3.4  A  [77],  respectively, 
we  expect  only  minor  differences  between  the  length 
scales  of  the  CO2-PH2  and  C02-He  potentials. 
Subtracting  the  1.9  A  effective  radius  of  the  pH2 
molecules  from  the  angle  dependent  minimum  of  the 
C02-He  potential  produces  the  outline  of  the  CO2 
molecule.  The  "dip"  in  towards  the  central  C  atom 
indicates  the  presence  of  an  attractive  well,  and 
correctly  reflects  the  experimentally  determined  T- 
shaped  global  minimum  C02-He  geometry  [78]. 
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Judging  from  the  drawings  in  figure  4,  we 
consider  it  unlikely,  but  not  impossible,  for  a  CO2 
molecule  to  occupy  a  single-substitutional  vacancy  in 
solid  pH2.  A  relaxed  trapping  site  structure  would 
require  considerable  distortion  of  the  surrounding  pH2 
molecules  away  from  their  original  lattice  positions, 
leading  to  numerous  smaller  repulsive  PH2-PH2 
interactions.  On  the  other  hand,  the  double- 
substitutional  site  readily  accommodates  a  CO2 
molecule,  and  the  triple-substitutional  site  appears 
perhaps  overly  spacious. 

We  can  most  easily  explain  our  observation  of 
three  main  peaks  in  the  spectra  of  as-deposited  CO2/PH2 
solids  by  trapping  of  CO2  molecules  in  three  distinct 
double-substitutional  vacancies.  In  each  case  the  CO2 
molecules  should  exist  as  strongly  hindered  librators 
and  not  as  rotors.  By  comparison  with  the  recent 
CO2/N2  experiments  [57],  we  would  expect  a  single 
sharp  line  as  the  main  signature  for  each  trapping  site. 
There  is  only  one  type  of  double-substitutional  site  in 
an  fee  solid,  while  there  are  two  different  ways  of 
removing  two  adjacent  host  molecules  in  an  hep  solid. 
These  are  designated  "in-plane”  (ip)  if  both  host 
molecules  are  removed  from  the  same  basal  plane,  and 
"out-of-plane"  (oop)  otherwise  [79]. 

Based  on  their  annealing  behaviors,  we  assign 
the  three  major  lines  to  CO2  molecules  trapped  in 
double-substitutional  sites  as  follows:  the  a  line  to  the 
oop-hep  site,  the  (3  line  to  the  ip-hep  site,  and  the  y  line 
to  the  fee  site.  The  virtual  disappearance  of  the  y  line 
upon  annealing  matches  our  previous  experience  with 
irreversible  loss  of  fee  regions  during  annealing.  We 
attribute  the  a  line  decrease  and  p  line  increase  upon 
annealing  to  the  improved  alignment  of  the  hep 
crystallite  c-axis  directions  with  the  substrate  surface 
normal.  Since  our  IR  beam  propagates  parallel  to  this 
same  direction,  CO2  molecules  oriented  perpendicular 
to  the  surface  normal  will  exhibit  better  alignment  of 
their  V3  mode  transition  dipoles  with  the  electric  field  of 
the  IR  beam. 

We  exclude  trapping  of  CO2  molecules  in 
triple-substitutional  sites  because  such  molecules 
should  exist  as  slightly  hindered  rotors.  Since  the 
rotational  constant  of  CO2  is  only  «  0.4  cm'\  several 
rotational  levels  would  be  populated  at  the  temperatures 
employed  in  these  experiments.  We  see  no  evidence  in 
the  spectra  for  the  evenly  spaced  absorption  lines 
expected  for  rotational  transitions  from  such  levels.  We 
note  that  our  results  do  not  absolutely  exclude  the 
possibility  of  trapping  of  CO2  molecules  ifrsingle- 
substitutional  sites.  In  analogy  to  the  argument  made 
previously  for  the  double-substitutional  sites,  one  can 
conceive  of  two  types  of  relaxed  single-substitutional 


sites  in  hep  pH2  which  differ  in  the  orientation  of  the 
trapped  CO2  molecule  with  respect  to  the  hep  c-axis. 

One  final  point  of  comparison  is  the  absorption 
linewidths  expected  for  the  single-  vs,  double- 
substitutional  trapping  sites.  In  the  first  case,  both  ends 
of  the  CO2  molecule  would  experience  strong  repulsive 
interactions  with  the  nearest-neighbor  pH2  molecules, 
suggesting  the  likelihood  of  strong  coupling  between 
the  CO2  V3  asymmetrical  stretching  mode  and  the  pH2 
solid.  For  the  double-substitutional  site,  especially  if 
the  center  of  mass  of  the  CO2  molecule  is  strongly 
localized  near  the  center  of  the  trapping  site,  the  ends  of 
the  CO2  molecule  should  experience  much  weaker 
attractive  interactions  with  those  pH2  nearest-neighbors 
aligned  with  the  molecular  axis.  In  this  case,  we  would 
expect  very  weak  coupling  of  the  V3  mode  with  the  pH2 
solid,  and  correspondingly  narrow  absorption  lines. 

C.  Weak  Absorption  Features 

We  turn  finally  to  the  origins  of  the  numerous 
weaker  absorption  features  in  the  CO2/PH2  spectra. 
Ultimately,  successful  quantitative  modeling  of  these 
features  should  yield  the  best  insights  into  the  CO2 
trapping  site  structures  and  librational  dynamics,  as 
well  as  the  more  general  lessons  we  hope  will  help  us 
interpret  the  spectra  of  larger  non-rotating  dopant 
systems.  At  present,  however,  we  can  only  offer  our 
qualitative  speculations  in  the  form  of  a  list  of 
hypotheses  to  be  examined  in  the  future. 

The  same  multiple  trapping  site  argument  used 
above  to  explain  the  a,  P,  and  y  peaks  can  be  extended 
to  include  some  or  all  of  the  peaks  observed  in  the 
2343.5  to  2345  cm’^  region  in  as-deposited  CO2/PH2 
samples.  In  place  of  double-substitutional  sites  in  hep 
and  fee  regions  we  substitute  the  concepts  of  hep- like 
and  fcc-like  trapping  sites  discussed  above.  In  contrast 
to  the  CH4/PH2  case,  for  which  we  found  no  such 
distinction,  for  CO2/PH2  subtle  differences  in  final 
relaxed  trapping  site  structures  due  to  non-nearest- 
neighbor  interactions  might  induce  the  ~  10  % 
differences  in  gas-to-matrix  shifts  required  to  account 
for  these  minor  peaks.  The  irreversible  changes 
observed  upon  annealing  would  thus  correspond  to  the 
disappearance  of  some  of  these  metastable  trapping  site 
structures.  The  logical  extreme  of  this  model,  in  which 
each  and  every  line  corresponds  to  the  absorption  of 
CO2  molecules  in  a  particular  trapping  site,  raises  the 
interesting  notion  of  inhomogeneities  as  discrete, 
countable  entities;  challenging  the  traditional 
"continuous"  depiction  of  inhomogeneous  broadening. 

The  extreme  dilution  of  our  samples  weighs 
against  attributing  the  weak  features  in  the  2343.5  to 
2345  cm’^  region  to  long-range  interactions  between 
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CO2  molecules.  Vibrational  shifts  due  to  resonant 
interactions  between  CO2  molecules  die  off  with 
increasing  separation  as  R'^  [46];  their  absolute 
magnitudes  fall  below  ~  0.1  cm"’  for  separations  greater 
than  w  20  A.  At  a  concentration  of  1  PPM  in  solid  pH2, 
the  number  density  of  CO2  molecules  is  only  2.6x1 0'^ 
cm'^;  for  a  random  distribution  the  mean  separation 
between  closest  CO2  molecules  is  «  200  A. 

Furthermore,  denoting  the  distribution  of  separations 
between  closest  CO2  molecules  as  f(R)dR,  then  from 
the  peak  value  near  R  =  200  A,  f(R)  tends  towards  zero 
as  R^  for  decreasing  separations.  Thus,  very  few  of  the 
CO2  molecules  will  be  close  enough  to  communicate  in 
this  manner.  An  intriguing  mechanism  in  which 
dopants  interact  via  the  strain  fields  they  induce  in  the 
surrounding  host  lattice  has  been  proposed  to  explain 
sharp  features  in  IR  absorption  spectra  of  SF6/Ar 
matrices  [80].  In  fact,  the  spectra  reported  in  that  study 
bear  a  passing  resemblance  to  our  CO2/PH2  spectra. 
However,  we  find  it  difficult  to  accept  that  such  a 
mechanism  could  operate  over  the  ~  100  matrix  host 
diameters  intervening  between  CO2  dopants  in  our  1 
PPM  samples. 

The  models  proposed  so  far  in  this  section 
cannot  explain  the  observed  reversible  temperature 
dependent  changes  to  the  spectra.  The  a"  line  in 
particular,  along  with  several  other  weaker  features, 
appear  to  originate  from  thermally  populated  initial 
states.  We  estimate  that  the  height  of  the  a"  peak 
grows  «  50-fold  upon  warming  from  T  =  2.4  to  4.8  K, 
implying  a  thermal  excitation  energy  of  ~  10  cm’’. 

Since  the  lowest  lying  vibrational  state  of  a  gas-phase 
CO2  molecule  is  the  V2  bending  mode  at  667  cm  ',  these 
thermally  populated  states  must  involve  interactions 
with  the  pH2  solid,  although  we  can  as  yet  offer  no 
specifics  as  to  the  libration-translation  dynamics 
involved. 

Another  possibility  is  the  aggregation  of  0H2 
impurities  around  the  CO2  molecules.  Assuming 
complete  equilibration  in  our  ortho/para  converter  [17] 
at  T  =  15  K,  the  residual  0H2  concentration  is  «  100 
PPM.  Actually,  assuming  natural  isotopic  abundance, 
the  pH2  solids  should  also  contain  »  300  PPM  HD 
molecules,  but  HD  molecules  are  immobilized  in  the 
solid  at  IHe  temperatures.  On  the  other  hand,  the  J=1 
"ortho"  excitations  are  mobile  by  a  process  known  as 
"quantum  diffusion,"  i.e.:  the  conversion  of  adjacent 
0H2-PH2  pairs  into  PH2-0H2  pairs  [69,73,81,82].  Any 
J=1  excitation  that  diffuses  into  the  vicinity  of  a  trapped 
CO2  molecule  would  localize  to  form  an  0H2-CO2  pair 
bound  together  via  the  electric  quadrupole-quadrupole 
(EQQ)  interaction.  Based  on  previously  published 
HD/0H2/PH2  and  D2/0H2/PH2  spectra  [4,5],  this  same 


EQQ  interaction  should  be  capable  of  producing 
complex  features  over  a  ~  1  cm'*  region  of  the  CO2 
spectrum.  In  this  picture,  the  reversible  temperature 
dependent  spectral  changes  could  be  attributed  to 
absorptions  involving  different  "orientational"  states  of 
the  0H2  molecules,  rather  than  the  libration-translation 
dynamics  of  the  CO2  molecules. 

Given  the  present  data,  we  can  only  note  in 
passing  that  the  weak  features  in  the  2345  to  2348  cm  * 
region  shown  in  figure  3  which  grow  in  upon  annealing 
show  blue-shifts  relative  to  the  a  and  p  peaks  consistent 
with  clustering  of  a  small  fraction  of  the  CO2  molecules 
to  form  CO2  dimers  and  trimers. 

V.  Conclusions  and  Future  Directions 

We  report  high  resolution  IR  absorption  spectra 
of  CO2  molecules  isolated  in  rapid  vapor  deposited  pH2 
solids.  The  absence  of  a  regularly  spaced  progression 
of  lines  is  taken  as  evidence  that  the  CO2  molecules 
cannot  rotate,  and  instead  exist  as  strongly  hindered 
librators.  The  three  main  absorption  features  are 
assigned  to  CO2  molecules  trapped  in  double- 
substitutional  sites  in  hep  and  fee  regions  of  the  as- 
deposited  solids.  Numerous  weaker  spectral  features 
remain  unexplained.  Reversible  temperature  dependent 
intensities  demonstrate  that  some  of  these  weaker 
transitions  originate  in  thermally  populated  states  lying 
~  10  cm''  above  the  respective  ground  states.  Some  of 
the  2  to  5  cm'*  blue-shifted  weak  features  observed 
upon  annealing  of  the  samples  may  be  due  to  small 
(C02)n  clusters. 

We  are  preparing  to  perform  a  new  series  of 
experiments  on  CO2/PH2  samples,  including  the  use  of 
isotopically  substituted  CO2  molecules.  Breaking  the 
inversion  symmetry  of  the  CO2  molecule  by  trapping 
the  ‘®0*^C**0  isotopomer  will  alter  the  nature  of  the 
allowed  librational  states.  Trapping  gradually 
increasing  concentrations  of  CO2  molecules  should 
permit  the  identification  of  spectral  features  due  to 
dimers  and  larger  clusters.  Manipulating  the  residual 
0H2  concentration  by  varying  the  ortho/para  converter 
temperature  should  clarify  the  importance  of  0H2-CO2 
clustering. 

We  expect  that  quantum  simulation  techniques 
will  be  very  valuable  in  elucidating  the  trapping  site 
structures  and  librational  dynamics  of  the  CO2/PH2 
system,  as  demonstrated  previously  for  O2  doped 
hydrogen  solids  [83-86]. 
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